NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 





1 v " : ■; 

(N&3A-CH-163 S/b6) UN THK itELATlONSHIi? N81-17177 

BETWEEN ENUINEERINU PHOPEETIES AND 
DELAMiNATiON UF COMPOBiTi: x^AtliNATES Interim 
Eeport (Vicyinid Polytechnic Inst, and Etate Oncias 


univ.) 25 P HC AU2/Mi- Aul CSUL llU GA/^'<4 44126 



COMPOSITE 


MATERIALS 


RESEARCH 


AND 


EDUCATION 


March 20, 1981 


Notice! Table 2 was inadvertently left out of the VPI-E-81-2 report you 
recently received entitled "On the Relationship Betxyeen Engineering 
Properties and Delamination of Composite Laminates." The table is 
etiiclosed for completeness. Please insert it after Table 1. 


THE NASA-VIRGINIA TECH COMPOSITES PROGRAM 

Nalional Aaronaulios & Space Adminisiration Virginia Polylochnic Insliluto 

Langley Research Comer and Stale University 

Hampton. Virginia 23665 Phone 703-961 -53 72 


reply to The NASA-Virginia Tech Composites Program 
Dept ol Engineering Science & Mechanics 
Virginia Tech 

Blacksburg. Virginia 24061 


N81-17177 


ON THE RELATIONSHIP BETWEEN ENGINEERING 
PROPERTIES AND DELAMINATION OF COMPOSITE 
LAMINATES 


C T. Herakovich 


Virginia Polytechnic Institute and State University 
Blacksburg, Virginia 


February 1981 


ON THE RELATIONSHIP BETWEEN ENGINEERING 
PROPERTIES AND DELAMINATION OF COMPOSITE LAMINATES 


by 

CARL T. HERAKOVICH* 


ABSTRACT : 

The influence of the coefficient of mutual influence, poisson's 
ratio and coefficients of thermal and moisture expansion on delamination 
is studied. Engineering theories are compared to finite element and expe- 
rimental results. It is shown that the mismatch in coefficient of mutual 
influence can have a strong influence on delamination with fiber angles 
in the 10°-15° range being critical for adjacent (±6) layer combinations. 
The mismatch in coefficient of mutual influence is reduced by a factor of 
two and the interlaminar shear stress is reduced significantly when 
the ±0 layers are interspersed between 0° and 90° layers. It is shown how 
the results can be used for design of composite laminates. 
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INTRODUCTION 


The interlaminar stresses in the boundary layer along the free edges 
of a laininated composite play an important role in the initiation of da- 
mage at the free edge. Delamination between layers is a direct result of 
these interlaminar stresses. While delamination is important to the ove- 
rall structural performance of laminates under tensile loading, it is 
critical for compression and shear loading where stability is a major 
concern. 

Following the initial work of Pipes and Pagano [ i ] , numerous inves- 
tigators have studied the problem of free edge effects in laminated com- 
posites. However, a systematic discussion of the fundamental relation- 
ships between the engineering properties of the individual layers and 
edge effects has apparently never been presented. This paper will attempt 
to at least partially fill this gap in the literature. 
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BASIC CONSIDERATIONS 


Mechanical loading : 


The fundamental reason for the presence of interlaminar stresses in 
laminated composites is the existence of a mismatch in engineering proper- 
ties between layers. For the purpose of discussion we consider here the 
special case of a balanced, symmetric coupon under axial loading in the x 
direction (Fig. 1) . The two most important lamina properties this 
problem are poisson's ratio 

~^Y “®i2 "" ^ ^12 ^ ^ ®22 ” 

( 1 ) V = -X = _ — = J= d , 

^ Sjj Sjj + (2 S,, + Sgg) + n“ 

and the coefficient of mutual influence 


( 2 ) n, 


'xy 




,X S„ 


'n m ^(2 Sjj ~2 Sjj-Sjg)^ in (2 Sjj ,-2 8 ^ 2 + S^g) 

.1 - I I «i y . n il.— .■■ M l. — II ,||.».W..||||.| 

j~m‘' Sjj + (2 Sjg + Sgg) + n“ SjgJ 


where S^^ and Sj^j correspond to the compliance coefficients in material 
principal and global coordinates, respectively, m and n are the oos9 
and sin0, respectively, and 0 is the angle of fiber orientation measured 
from the axis of the coupon. Equations ( 1 ) and ( 2 ) are plotted in Fig. 2 
for properties typical of T300/5208 graphite/epoxy. It is noted that 

is an odd function of 0 with a maximum value at 0 = 11,5“ whereas 


'xy 


X 


V is even in 0 with a maximum value at 0 = 22°, 
xy 


If there is no mismatch of v or ri,„ between layers, there are no 

^ ^,x 

interlaminar stresses regardless of the mismatch in elastic and shear mo- 
duli. The mismatch of poisson's ratio between adjacent layers gives rise 
to dissimilar lateral strains E^ in free (unbonded) layers, but results 
in identical strains at the layer interfaces with accompanying interlami- 
nar stresses a and t in perfectly bonded laminates. Likewise, the 

y 
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ORIGINAL PAGH US 
OF POOR QUALITY 

n mismatch gives rise to nonzero interlaminar shear stresses T in 
Xy ^ zx 

bonded laminates. The magnitude of the interlaminar stresses is related 


to the magnitude of the mismatch in V and r|, 




the elastic and shear 


moduli and the stacking sequence. The stacking sequence plays an impor- 
tant role in that is establishes the moment associated with the stres- 
ses. This "loment is equilibrated by the interlaminar a distribution 

z 

(Fig. 1). 


In order to demonstrate the fundamental relationship between engi- 
neering properties and the interlaminar stresses, we consider three spe- 
cial classes of adjacent laminae ; { +0/-0 ) , (0/0 ) , and ( 0/90) . The 
properties of greatest interest are the mismatch between adjacent layers, 


6v and 

xy 

with 


( 3 ) 

1 6v 1 = 

1 V (0„) - V (01 ) 1 

' xy xy ^ ' 

( 4 ) 

1 1 
^y,x 

= 1 n (9,) - n (e,) 

' xy 2 'xy 1 

, X , X 


where 0j and 0g correspond to the fiber orientations of adjacent layers. 
Interlaminar stresses will be largest when the mismatch is maximum (all 
other things being equal) . 


Fig. 3 shows plots of | <Sv | and 1 6ri„, 

xy xy 


as a function of 0 for the 


X, 


three special cases under consideration. Thd largest mismatch is observed 


'xy 


= 4.34 for the (+ 11.5/-11.5) combination. The curve exhi- 


X 


to be I (Sn 

bists a ste^p gradient in the vicinity of the critical value 0 = 11.5°, 

The magnitude of <Sri for the (+0/-0) combination is exactly double that 

X 

for the (0/0) and (0/§O) combinations. This is because n , is an odd 

^.x 

function of 0 and n = 0 for 0=0° and 0 = 90° . There is no mismatch 

^,x 

in V for the (+0/-0) combination as V is an even function of 0. 
xy *^0; 


The majcimum value of l<5v^| is much smaller than 1 for all 

three laminae combinations. The worse case for 6v is the '^0/90) combi- 

xy 


nation with 6v 


xy' 


0.34 at 0 = 22°. Thus, the maximum value of 6n. 


is more than ten times the maximum 


, I ^ X 

[(Sv |. Since the interlaminar normal 
xy 
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and shear strengths are of the same order of magnitude, initiation of de- 
lamination is much more sensitive to the mismatch in coefficients of mu- 
tual influence than the mismatch is poisson's ratio, (as will be shown 
later, this large difference in mismatch can be partially or totally of- 
fset by the influence of stacking sequence on the magnitude of a ) , 

z 

close examination of Pig. 3 shows that |<So^y is significantly 

larger thanliv) | for almost the entire range of fiber angles for all 
xy 

three laminae combinations, the exception being as the (0/90) combination 

(with “ 0) is approached. The zero value of and the near 

maximum value of for the (0/90) combination is consistent with the 

rather well-know numerical result that T is zero in crossply laminates 

zx 

and delamination is solely due to interlaminar normal stresses. 
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OF POOR QUAtrfV 


Thermal and Moisture Effects 


The mismatch of the coefficients of thermal and hygroscopic expan- 
sion produce edge effects similar to those due to mechanical loading [2], 
The formulation of both problems is identical and hence we consider only 
the thermal problem for the purpose of illustration. The free thermal 
expansions for a uniform, unit temperature change of an unbonded lamina 
are 

rp 0 •> 

I \ I rrr + n a.2 \ 

( 5 ) < ? = I a, + a, [ 

^ y I i » 

\ / \ 2 m n (ai - 02 ) / 

and the mismatch in strains of interest at the free edge can be expressed 

( 6 ) 6 = e!^,( 2) - ej(l) = (ag - a^) (cos^ 0z - oos^ 0i) 

y y y 

~ “1^ 2 0^ - sin 2 02) 

xy xy Ay 

In the above, a's are coefficient of thermal expansion, subscripts 
X and y refer to global coordinates, and the superscript T indicates ther- 
mal effects. 

Application of eqns. (6.7) f.or the three laminae combinations under 
consideration gives the following results : 

(±0) Combination : 

fiEy = o 

( 8 ) 6y^ = (ttj “ aj) (2 sin 2 0) 

|6y^| = 2 (02 - ai) at 0 = 45° . 

max 
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(0/0) combination 


OF POOR QUALfiT 


fiEy = (ttj - aj) (sin^ 0) 

|6e^| = (ttj - otj) at 0 = 90“ 

^ max 

(Sy'^, = (aj - tti) (sin 2 9) 

= (a 2 - aj) at 0 = 45® . 

^ max 


(0/90) Combination : 

SSy = (a2 - aj) (“Gos^ 9) 

I'^^ylmax = ^""2 " 

I 'Sy'LI = (aa - ai) (sin 2 0) 

xy 

I<Sy^I = (aa - ai) at 0 = 45° . 

^ max 


T 

It is evident from the above that the worse case is 6y for the 

T ^ 

(±0) combination with 0 = 45°. The magnitude of Sy^,, for the (±0) combi- 
nation is twice that of the others, as it was for 6n j however, the 

^,X ip 

critical angle has shifted from 11.5° to 45°. The mismatch is zero 

for the (-0) combination, as it was for 6v , but the maximum mismatch 
for the other two layer configurations both correspond to the (0/90) case 

unlike 6v where the maximum mismatch occurred at 0 = 22°. It is also 

T T 

noted that the maximum 6y is only double that of the maximum 6e com- 

'xy ^ y 

pared to a factor of more than ten between 1 6n | and | <Sv | 

iiTax ^ max 

Also, the variations in thermal strain mismatch as a function of 0 are 
smooth trigonometric functions and do not exhibit the sharp rise to a peak 


value as does 6f| 


'xy 


in Fig. 3, 


,x 
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Since composites are cured at an elevated temperature thermal stres- 
ses are always present. For a negative change in temperature (as during 

T 

cure) the thermal shear strains have the opposite sign as those due 
to tensile loading and hence the effects tend to offset, as well as shift 
the critical angle. The thermal and mechanical effects would, of course, 
be additive for compressive loading. The combination of thermal and mecha- 
nical loading effects on iSe^ is additive for tensile loading and offset- 
ting for compressive loading. 


NUMERICAL RESULTS FOR I ± 0 1 5. LAMINATES 


The discussion in the previous section is based upon somewhat heuris- 
tic arguments. Fig. 4 shows finite element results for the influence of 
fiber orientation on the value of the interlaminar shear stress at the 

intersection of the interface and the free edge in [±0] laminates [3l. 

s 

The nondimensionalized stresses in the figure correspond to the maximum 
stress for each laminate configuration. The results were obtained for selec- 
ted fiber angles of T300/52O8 graph! te/epoxy. The curve clearly shows a peak 
value at 0 = 15°, close to the value of 11.5° for the maximum |6n. 


and there is a distinct similarity between the distribution of [iSri 


xy,x' 


j<y 


in 


Fig. 3 and T in Fig. 4. 


X 


A failure analysis for the [±0] laminates was performed using the 
finite element stress results and the Tsai-Wu tensor polynomial failure 
criterion [ 4 ] . Figure 5 shows the variation of the maximum value of the 
tensor polynomial in these laminates as a function of a fiber orientation. 
Again, there is a distinct peak value at 0 = 15°. Examination of the indi- 
vidual terms of the polynomial clearly indicated a dominate influence of 
T for this critical fiber orientation. Thus, there is numerical support 

2X 

for the heuristic arguments of the previous section. 
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EXPERIMENTAL RESULTS FOR [ - 0 Ig LAMINATES 


Tho mode of failure at the free edge of [10]_ Colion 6000/PMR-15 

s 

graphite/epoxy laminates was determined in reference [ 5 ] using the edge 

replicating technique [ 6 ] . Figs 6-8 show replicas of the edge damage 

for selected laminates. The mode of damage is delamination for [ilo]_ 

s 

and [ ±30 J „ laminates, but shifts to transverse cracking for the [±60]_ 
s » 

laminate. These results give experimental support for the dominate influ- 
ence of ti in the 10°-30° range of fiber orientations for [±0]„ lomi- 
nates. 


LAMINATE DESIGN 


The results in Fig. 3 can serve as a guide to the designer in choo- 
sing fiber orientations and stacking sequence in applications with free 
edges. Consider the case of the quasi-isotropic laminate with fiber orien- 
tations of 0°, 90°, + 45° and -45°. An often asked question is : what is 
the optimum stacking sequence ? The curves in Fig. 3 suggest that the in- 
terlaminar shear stresses T will be significantly lower if the + 45° and 
-45° layers are separated by a 0° or 90° layer. 

Of the twenty-four possible combinations of layers in this laminate, 
only twelve are distinct because of the interchangeability of the +45° 
and - 45° layers. Of these twelve, there are six possibilities with adja- 
cent ± 45° layers and six with ± 45° layers interspersed between 0° and 
90” layers. These twelve laminates are depicted in Tables 1 & 2 along with 
the magnitude and direction of the 0^ stresses for tensile loading £^^=0,5%. 
Also indicated in the Tables is pertinent information as to the magnitude 
and direction of the equilibrating moment at each interface, and the maxi- 
mum a and T stresses determined from a finite element stress analysis. 

z zx 
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A hierarchy of lamiantes for resistance to delamination can be esta- 
blished through consideration of the magnitude of the magnitude 

and direction of the interface moment. The hierarchy given in the Tables 
is based upon the following guidelines listed in order of senverity : 

1. avoid adjacent ± 45° layers 

2. select the most negative or smallest positive interface 
moment 

The analysis indicates that the [ 90/45/0-45 ] laminate provides the 

s 

greatest resistance to delamination for tensile loading. This laminate 
has the ± 45° layers interspersed and the interface moment is negative 
throughout. In contrast, the [ 45/0/-45/90 ] laminate is the worst case 
of the six with interspersed ± 45° layers, because the interface moment 
is always positive and of maximum magnitude. 

The hierarchy presented in the tables does not include the influence 
of combined stress states nor the effects of residual thermal or hygros- 
sopic stresses. Thus it is expected that experimental results may indicate 
some reordering of the list. A linear elastic finite element stress analy- 
sis of the twelve possible laminate configurations was conducted to obtain 
approximate elastic stress distributions. The magnitude and location of 

the maximum values of a and T are presented in the tables along with a 

z zx 

strength hierarchy based upon the following guidelines : 

1 . minimize I t I 

' zx' 

2 . minimize a 

z 

The finite element results clearly show the advantages to be gained 
with interspersed ± 45° layers. All six laminates with interspersed ± 45° 
layers have lower interlaminar shear stresses than those with adjacent 
* 45° layers. The magnitude of the interlaminar shear stress ranges from a low 
of 5.8 KSI (39.9 MPa) for the [ 45/O/90/-45 ] ^ laminate to a high of 
9.2 KSI (63.8 MPa) for the [ 90/45/-45/0 ] laminate , a sixty percent dif- 

O 

ference. The range of a stresses for interspersed ± 45° layers is smal- 

z 

ler than that for the adjacent layer configuration indicating a more ef- 
ficient design with interspersed ± 45' sin the presence of alternating 
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positive/negative loads. The hierarchy based upon the finite element stress j 

results is the same as that based upon the engineering or heuristic ap- j 

proach . | 

i 

It is interesting to note that the great majority of experimental | 

results reported in the literature for delamination studies of quasi- i 

isotropic laminates have been obtained from specimens with adjacent ± 45° 
layers [6-8]. The only exception known to the author is that of reference 
9 where it was shown that the strain at first delamination in a j 

[145/0/90] laminate was approximately half the strain at first delamina- i 

tion in a [ 0/45/90/-45 ] laminate. The [1 45/0/90 ]_ specimen delaminated i 

S S' I 

along the midplane due to high interlaminar normal stresses and the | 

[ 0/45/90/-45 ] specimen delaminated in the 90° layer, apparently due 1 

to combined a„ and effects. I 

2 zx 

Further evidence of the importance of adjacent ± 45° layers on failure ; 

of quasi -isotropic laminates were presented in reference [6]. The final ! 

fracture patterns [ 0/145/90 ] ^ and [ 0/90/145 ]„ graphite-epoxy laminates 

s s 

exhibited extensive delamination along the 1 45° interface of both lamina- 
tes even though the interlaminar normal stress is compressive along this 
interface in the [ 0/90/145 ]^ , The finite element results of this study in- 
dicated very high interlaminar shear stresses for this interface (Table 2) . 
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DISCUSSION 


Do lamination has been studied with an emphasis on straight coupons 
under tensile loading. It has been pointed out that adjacent * Q layers 
leads to high interlaminar shear stresses and that interlaminar normal 
stresses are very dependent on stacking soquoneo. Since shear failure io 
independent of loading direction, adjacent i 0 layers should generally 
bo avoided. There are a limited number of situations whore the free edges 
are al^^;ays parallel to the loading axis and the direction of loading is 
know to be constant. In such situations tliore may bo justification for 
adjacent f 0 layers in order to optimisso the negative interlauninar normal 
stress ; however this can only bo accomplished at the expense of signifi- 
cantly higher interlaminar shear stress (Table 2) . In many application the 
loading direction will change during service and free edges are often pre- 
sent in the form of holes and cutouts eliminating the . ondition of a single 
known direction of loading. Interspersed f 0 layers are preferred in such 
applications as inter laiininar shear stresses are reduced as are the extre- 
mes of interlaminar normal stresses. An example has boon given as to how 
these results can bo used systematically in design to rainimiae the pos- 
sibility of dolamination. 
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CONCLUSIONS 


It has boon shown that thoro As a fundamental relationship between 
dolnmination and the "mismateh" of engineering properties between adja- 
cent layers of laminated composites. Interlaminar shear stresses are 
primarily a function of the mismateh in coefficients of mutual influence 
which can bo as much as ton times greater than the mismatch in poisson's 
ratio. The mismatch in coefficients of mutual influence has a high peak 
value in the fiber range for f 0 laminae combinations. This mis- 

miatch is reduced by a factor of two when the ^ 0 layers are interspersed 
between 0° and 90“ layers. Interspersed i 0 are recommended for more ef- 
ficient design. 
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Figure 7 : EDGE DAMAGE IN A [(±3i CELION 6000/PMR-15 

L S 

GRAPHITE-POLYIMIDE LAMINATE AT 99% OF ULTIMATE STRAIN 
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Figure 8 : EDGE DAMAGE IN A [(± 60 ) 2 ]^ CELION 6000/PMR-15 

GRAPHITE-POLYIMIDE LAMINATE AT 99% OF ULTIMATE STRAIN 
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Table 1 ; Strength hierarchy for interspersed ±45“ layers. 























o Maximum T300/5208 Graphite/Epoxy 

D Maximum lo I E = 0.5% 

' z ' X 

Table 2 : Strength hierarchy for adjacent ±45° layers. 

























